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ABSTRACT 


The primary object of this report is to develop a linear model that 
describes the perturbation motion of a spinning spacecraft made of a family of 
tanks that are partially filled with fluid. 

First the fundamental nonlinear equations of motion are derived and then 
specialized to a steady-state rotation of the vehicle about a given axis of 
rotation. Then, a thrust about the spin axis is introduced. Finally, a pertur- 
bation solution is derived which linearizes the problem. The effect of the 
centrifugal and corlolis accelerations together with vorticity are Implicitly 
taken into consideration in the formulation. A variational formulation of the 
associated boundary conditions is presented. For most practical cases it is 
shown that the simple classical pendulum representation for slosh is not very 
appealing for a spinning spacecraft unless severe restrictions are allowed. 
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INTRODUCTION 


Problems concerning the behavior of liquid propellant motions In 
partially filled tanks of spin-stabilized spacecraft have gained topical 
importance recently, especially In those Instances where the propellants con- 
stitute a large percentage of the total vehicle mass. As is known, the 
Interaction between contained propellant and rigid body motions plays a sig- 
nificant role In establishing the design of the attitude control system and 
assessing the performance of a spacecraft. Indeed, the action of liquid pro- 
pellant motions on the tank walls ("sloshing") of a spin-stabilized space- 
craft may couple with Its natural nutational mode of motion so as to saturate 
the control system and possibly cause flight Instabilities, or severely 
impair Its performance such as maintaining pointing accuracy. To predict and 
control the mot.ons of a spin-stabilized spacecraft with partially filled 
liquid cavities, it is therefore necessary to single out the participation of 
sloshing propellants in the overall vehicle motions. Methods for character- 
izing the behavior of propellant motions In partially filled spherical tanks, 
as encountered With the design of the Galileo Dual-Spin .Jupiter Orbiter space- 
craft, tfore originally developed under somewhat restrictive assumptions. 

In 1978 Ref. 1 devised a perturbational construct describing the small 
liquid motions in a partly filled spherical tank relative to a known refer- 
ence state of motion. The reference motion was defined as one of uniform 
rotation of the liquid cavity about a fixed point on a stationary axis under 
the action of a constant acceleration directed negatively along that axis. In 
arriving at this construct, it was assumed that: 1) the propellant was an 

Inviscld fluid; 2) the vortlclty was Independent of the spatial coordinates, 
i.e., a function of time only; 3) the steady state spin and all rotational 
perturbations were instantaneously communicated to the liquid; 4) Coriolis 
accelerations were significant throughout the liquid volume but negligible 
at the free surface; 5) translational motions of the liquid tank systems were 
nonexistent. 

In 1979 Ref. 2 showed that assumptions 2) and 4) gave rise to spurious 
instabilities for certain liquid fill ratios, and that assumption 3) was incon- 
sistent with the mechanism of vortex transfer from the tank wall to the 
liquid. The Inconsistency of the construct was further substantiated by 
Ref. 3. 


The primary aim of this investigation was to develop a linear construct 
that describes the perturbational motions of a spin-stabilized spacecraft, 
comprising N tanks partially filled with liquid propellants, relative to a 
given state of motion. A further objective was to establish conditions under 
which sloshing motions in spinning tanks could be represented by simple 
mechanical pendulums. 

The basic subject matter is kept as concise as possible. In Section 1, 
certain fundamental results and equations from classical mechanics and hydro- 
dynamics are used to derive the equations of motion of a spin-stabilized 
spacecraft comprising N tanks partly filled With liquid propellants. The par- 
ticipation of the motions of the liquids in the overall vehicle motion is 
singled out. Diverse forms of the equations are presented. The fundamental 
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equations are specialized to the steady-state rotation of the vehicle, as a 
single rigid body, about its designed spin axis under the action of a con- 
stant thrust directed along that axis. In Section 2, the nonlinear equations 
are specialized to a reference state of motion in which the spacecraft 
(solid + liquids) undergoes a uniform steady state rotation about its 
designed spin axis under the influence of a constant thrust directed along 
that axis. Then, the vehicle Is slightly disturbed from the known reference 
state of motion and the resulting perturbational equations of motion obtained. 
The perturbational equations are modified to reflect the correct mechanisms 
of vorticity transfer from the boundaries of the tanks to the contained 
liquids in accordance with Ref. 2. In Section 3, the perturbational liquid 
motions are made to depend on the solutions of N Inhomogeneous boundary value 
problems. The effects of centrifugal and coriolis accelerations together 
with vorticity are implicitly taken into consideration In the formulation. 
Variational formulations of the associated boundary value problems are pre- 
sented. For most practical applications, It is shown that the use of mechan- 
ical pendulums to represent sloshing motions in spinning tanks is suspect and 
should be avoided, unless severe restrictions can be tolerated. In Section 4 
the resulting perturbational equations are summarized. 
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SECTION 1 


EQUATIONS OF MOTION OF A SPIN-STABILIZED SPACECRAFT 
COMPRISING N TANKS PARTLY FILLED WITH 
UQUTlt PROPELLANTS 


PHYSICAL ABSTRACTION 

Imagine a spin-stabilized spacecraft carrying N tanks partly filled with 
liquid propellants, as shown in Figure 1* It Is assumed that: (1) the 

liquid cavities are arranged In a symmetrical pattern about the designed spin 
axis of the vehicle; (2) In the abr-ence of liquids, the mass center of the 
spacecraft lies on its designed spin axis; (3) the propellants are Ideal 
liquids; (4) capillary and mass forces acting on the liquid are negligible; 
and (5) mass forces acting on the solid are Insignificant. 


COORDINATE SYSTEMS 

To describe the motion of the vehicle, take a righ£-hand£d cartesian 
frame of reference with origin <> and base vectors (T, j, k) fixed rela- 
tively to the solid, To complete the specification of assume that (’ Is 
located at the mass center of the spacecraft sans liquids, and that k Is 
coincident vlth Its designed spin axis. 

To describe the liquid motions in the nth tank, take a right-handejd 
.cartesian frame of reference /? n with origin C' n and base vectors (i n , j n , 
k n ) fixed relatively to the cavity, such thjat ,^n J| «iF. The position of ^'n 
relative to (> is designated by the vector L n , 

Finally, the motion of with respect to inertial space Is specified 
by the translational velocity vector u and angular velocity vector m. 


K1NEMATICAL MOTIONS 

Let r denote the position vector of a mass point In the solid relative 
to ijF , Then, the velocity and acceleration vectors of the point are simply. 


V * u + t\\ X r , 


= V + ojXV = a + wXr + m X (w X r) 


( 1 . 1 ) 


where 


a = u + w X u 


( 1 . 2 ) 


l-l 


stands for the acceleration of the spacecraft. It should be appreciated 
that differentiation with respect to ,‘F (and *^) Is Indicated by (.!.) 

Specify by r n the position vector of a material point Q n In the nth 
cavity relative to »^ n , l,e., relative to the cavity. Then, the velocity 
V n » V n (q n ,t) and acceleration ja n * a n (Q n ,t) of the liquid are given by 


V 


n 


where v n » 
i . e . , 


u + x R + « 
n 


n 


• V + ••• XV * 4 + r.' 
n n 

Vn (Qn.b) denotes the 


x R n + ~ x (“ x R n ) + 2” x v n + v n , 
velocity of the particle relative to 


(1.3) 




n» 



r 

n 


I 


and 


R 

n 




(1.4) 


(1.5) 


Note that Q n may be a point in either the liquid volume t n , free liquid sur- 
face nj? n or wetted boundary *i^ n of the tank, (See Figure 2 for notations.) 

The material derivative of a vector or scalar point function associated 
with the motion of liquid in the nth container is defined as 


(••*) * 3 ( ■ . • ) + (r 

.it 


n 


V ) (• 
n 


,) » (• 


,} t + 


(v_ 


V > (• 
n' 


where, of course 


V = i - 2 - + j — 2 — + k — L„ 2i„l, *.,‘J * (< 

n n flx 3 n By n Hs ’ 9t v 
n J n n 


• i • 


In particular, 


v=v + (v *V ) v = v + V 4 v^ + F, V, 
n n c n n n n t n 2 n n n 


( 1 . 6 ) 


(1.7) 


( 1 . 8 ) 
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whore £ n * Cn (Qn» fc ) denotes the vortieity tor of the liquid relative to 
**n» *•••» 


V x v 
n 


n 


(1.9) 


The vortlclty vector of the liquid In the nth tank relative to Inertial 
Space Cn * Cn (Qnit) is defined by 


'n 




2w + § 


n 


( 1 . 10 ) 


whlch_follows by taking the curl of (1.3i). For reference, note that If 
Vn * V n " V n * v n m 0 everywhere, then it can be shown that 


v x a * r, + « x c - (C • V ) V 
n n n n w n n' n 


'n - (5 n ’ V n > \ ~ 2 <“ ’ V v „ + 2 “ 


+ V „ X X V „) 

n. n n n 

t 


E + V X (? X v ) + 2w X f; - 2V (w • V ) + 2m 
n. n ’n n 'n n ' n 

t 


( 1 . 11 ) 


EQUATIONS OF MOTION 


Considering body and liquid as a single mechanical system subject 
to given forces, the equations of motion can, in view of the stated 
assumptions, be written 


N 


Ma » F + 




n=l 


N 

I • w + m x(I • w) * T + T t 

n=l 


0 . 12 ) 
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» 


— - V p 
P n*n 
n 


V * V • v 
n n n 


0 , y Q n e , (n*l,2, • ' • »N) 


V X a = 0 , 
n n 


(1.13) 


where (F, T) are the resultant vectors of external forces and moments 
applied to the system, (Fn» T n ) are the resultant vectors of forces and 
moments arising from the action of the liquid on the nth cavity wall, 
p n (Q n »t) and M n denote the pressure and density of the liquid in the nth 
tank, M stands for the mass of the solid, and I is the momental dyadic of 
the solid, l.e., 


I * f [R 2 E 
solid 


rr] dm , 


(E “ 11 + jj + kk) . 


(1.14) 


Clearly, equation (1.12) govern the translational and rotational motions of 
the vehicle. Equations (1.13), together with certain boundary conditions 
to be established, govern the motions of the liquids in the N cavities. 

For a particular tank, (1.13i) is the classical Eulerian equation of motion 
of a liquid in the absence of mass forces, and (1.132) I s c h e mass conserva- 
tion principle for an incompressible liquid, i.e., continuity equation. 
Another form of the mass conservation principle used implicitly throughout 
this Investigation is given by 


(p dt )’ = 0 , (n=l,2, »N) . (1.15) 

n n 

Equation ( 1 . 1 3 3) is the result of taking the curl of (l.l3i), l.e., the com- 
patibility (vortlclty) equation. 


BOUNDARY CONDITIONS 

The liquid velocity relative to each cavity must satisfy certain kine- 
matical conditions, namely, the component normal to the wetted surface a Wn 
of the tank must vanish (condition of non-adherence), and the component 
normal to the free boundary Op must equal the velocity of the surface 
normal to itself. Thus, if the equation of the free boundary is specified 
in the general form F n (Q n ,t) = 0, Q n e Op , one can write 

n 
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0 


> 



* 



Q e a 


Q e a„ 
^n F 


(n-1,2, • • • ,N) 


(1.16) 


where v n Is the unit vector along the outward drawn normal to the surface 
under consideration. 


The free boundary condition can also be written in the form 


F + v • V F = 0 
n. n n n 
t 


Q e a_, 
n F 

n 


(n=l,2, • • • ,N) 


(1.17) 


which follows by setting the material derivative of F n to zero. Clearly, 


*n- 7 nVI V n F J’ Q n e *F ’ (n=l , 2 , • • • ,N) . (1.18) 

n 


In particular, if 


? n = z n - f^ (x^, y n , t) =0, Q n e ct f , tn=l , 2, • ■ • ,N) . (1.19) 

n 


then the kinematic conditions at c/p become 

n 


v • v = f cos (v , k ) , Q e a , (n=l,2,..«, N) . (1.20) 

n n n_ n n n r 

t n 


Forms (1.19) and conditions (1.20) are used extensively later. 

Kinematical conditions (1.16), in conjunction with (1.132) and the 
divergence theorem yield the injunctions 



( 1 . 21 ) 


i.e., global continuity. 
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The pressure inside each cavity must satisfy a condition of thermo- 
dynamic equilibrium at the free boundary, which, in the absence of capillary 
forces, can be taken 


p n * 0 > Q n e cr p , (n=l,2, • • • ,N) 

N 


( 1 . 22 ) 


FORGES AND MOMENTS RESULTING FROM ACTION OF LIQUIDS 

The forces and moments (F n ,T n ) resulting from the action of the liquids 
on the walls of the cavities are given by 


r 


F =* lv p da = lv p do + I v p da„ 
n I n *n w J n r n w J n *n F 
^ n n ^ n 


w 


w 


n 


/ v p da = / V p dt , (n=l,2, • • • ,N) , 

n n n J n n n 


T * 
n 




f 


w 


( R t1 XJ) P n da 
n n n w 


/ 


w 


(R X v) p d 
n n w 


n 


n 


(1.23) 


+ / (R X v ) p 

J n n n 


da T 


n 


/ (R X v ) p da = (R X ? p ) dt , (n=l,2, • • • ,N) , 

n n *n n J n n r n / n ’ ’ ' 


n 
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in view of (1.22) and the divergence theorem, where 


0*0 + 0 _, , 
n w F 
n n 


(n*l,2, • • • ,N) 


(1.24) 


INSTANTANEOUS MASS CENTER OF SPACECRAFT 

Let the instantaneous position of the spacecraft mass center (solid + 

liquids) relative ro dF be specified by the vector r ; then 

c 


r » 

c M+M 


4r- f R p dx 

■Mj J n F n n 


(1.25) 


n=l x 


in light of 


/ 


r dm * 0 , 


solid 


(1.26) 


where the total mass of liquids 


N 


M, = 7 / p dx 

n n 

n=l x 

n 


(1.27) 


OTHER USEFUL FORMS OF THE EQUATIONS 

Equations (1.32), (1.13i), (1.23) and (1.25) make it possible to write 
(1.12) as 


N 

(M+M.) [a + u X r + wX(ujXr)] + ^ / (v + 2m X v )p dx = F , 

1 c c J n n 'n n ’ 


n=l x 


n 


(1.28) 


N 


K *ui + iiX(K *w) + (M+M x ) (r c Xot) + ^ J ( R n x (v n +2wX V n )]p n dx n = T, 

n=l x 

n 
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where the dyadic 


K - X + 


N 

E / 


n*l t 


n 


'l»„l B-W % dT n 


(1.29) 


Solving for a from (1.28^) and substituting the resulting expression Into 
(1 . 28 2 ) > 


N 

e + ;“ X ( 0 -r") + J [ (R n -r c ) X (v n + 2» X v n )]p n dt n - T - r, X F , 

n*l x 

n (1.30) 


where 

0 * K - (M + Mj) [|rj 2 E - r e r c ] . (1.31) 


Equation (1.30) clearly represents the rotational equation of motion of the 
vehicle about its instantaneous mass center. Note that I Is the momental 
dyadic of the solid whose components have constant values with respect to .P . 

K - I is the momental dyadic of the liquids whose components are functions 
of the instantaneous forms of the liquids inside the tanks. 0 -K is the 
inertia tensor whose components are functions of the Instantaneous position 
of the vehicle mass center relative to ,f. 

Diverse forms of (1 . 1 ) are also possible. In particular, 

v +2wXv +? xv + m X r + V «J> =0, (n=l , 2, • • • ,N) (1.32) 

n. n n n n n n 

t 


where the acceleration potential 


~ (w x R ) 2 , (n=l , 2, * • • ,N) , 

2 n 

(1.33) 


n 


r + H + Ca + i X V ‘ \ 

n 
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can be derived from (1.3)» (1.6) and (1.132). Taking the curl of (1.32) 
gives the compatibility equation 


F, + 2« X £ + V X (5, 


n 


n 


n 


n 


xv n ) 


2V (w 


v ) 
n 


+ 2w = 0 


(n“l , 2 N) . 

(1.34) 


STEADY STATE ROTATION OF VEHICLE UNDER ACTION OF CONSTANT THRUST 

If the applied forces and constraints are such that the spacecraft 
(solid + liquids) performs a uniform steady state rotation to as a single body 
about its designed spin axis k under the action of a constant thrust F 
directed along k, the system is In equilibrium relative to the moving frame 
of reference JF attached to the body, and (1.28), (1,30) can be written 


(M + M L ) [o + u X (5 X r c )] - F , 
to X ( 0 *<o) ■ T - r X F ** 0 


(1.35) 


But (1.352) Implies 


0 •{.' » C 


T - r X F - 0 
c 


(1,36) 


where c = c k. Equation (1.36i) states that the angular momentum about the 
center of mass is conserved. If 


r - r k , 
c c 


(1.37) 


then (1.362) Is satisfied Identically providing 


T S 0 . 


(1.38) 


With to = tok, r = r it and F = F k, equation (1.35,) yields 

O G •*’ 


M + M, 


k = sk , 


(1.39) 
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i 


The area parameter c in (1.36i) is given by 


c ■» 1 


kk 


(1.40) 


where 


I. , « k • 0 • k 
kk 


(1.41) 


is the moment of inertia of the spacecraft about the axis k. 
Now v^ * 0 in each tank, and (1,32) > (1.33), it follows 


1 n - - 1 - 2 

— + n • R - — (>,' x II ) *<j> ** c ,Q ct , (n*l ,2, • • • ,N) , 

p n 2 n n n ti n ' 

n 


(1.42) 


where c n are constants. If the equations of the free boundaries are speci- 
fied as F„ * z - f « 0, then (1.42) can be written 
11 n n 


p “ up (f - z ) , Out, 
r n n n n ui n 


(n*l ,2, • • • ,N) , 


(1.43) 


in light of conditions (1.22), where 


f n = a {c n + 1 ^ C(L ln + V 2 + (L 2n + y n )2]} + L 3n » 


(n=l ,2 , • • • ,N) , 
(1.44) 


and (L]_ n » 1.2n> l* 3 n ) are the components of L n with respect to tF. The 
equations 


F =z -f =0, Q e , (n=l,2, • • • ,N) , 

n n n hi F 

n 


(1.45) 


where f n are given by (1.44), are N paraboloids of revolution whose axes 
coincide with the vehicle spin axis. The constants c can be determined from 
the known liquid volumes. n 
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SECTION 2 

SMALL PERTURBATION EQUATIONS OF MOTION OF A SPIN-STABILIZED 
SPACECRAFT COMPRISING N TANKS PARTLY FILLED WITH LIQUID PROPELLANTS 


GENERAL CONSIDERATIONS 

The equations that govern the motion of a spacecraft carrying N tanks I 

partly filled with liquid propellants are nonlinear, and do not admit general, j 

exact solutions. Nevertheless, It is sometimes possible to obtain or approxi- { 

mate particular solutions, such as steady state rotation of the vehicle about its I 

designed spin axis under the action of a constant thrust derived above. Then, if | 

the system is supposed to be slightly disturbed from this known state of motion, | 

the resulting small motion of deviation can usually be described by a set of jj 

linear ordinary and partial differential equations with appropriate Initial and ;j 

boundary conditions. jj 

Such equations are usually constructed in the following manner. Suppose \\ 

the complete nonlinear equations of motion of the vehicle and liquids to be ■ 

written for two states of motion: reference and disturbed. The reference 

or 'undisturbed 1 state describes the path of the mass center, attitude of 
the vehicle, and the shapes of the free liquid surfaces in the cavities. The 
disturbed state describes the actual motions of the vehicle and liquids, and 

consists of the reference motions plus small deviant motions therefrom. The -j 

difference between the governing equations of the perturbed and reference j 

states, upon neglecting all terms involving products of deviant quantities i 

and/or their derivatives, yield': the perturbation equations (Ref. 4). 


REFERENCE STATE OF MOTION 

The reference state of motion embodying the particular desired attitude 
of the spacecraft is, in general, one of nonuniformity. For many practical 
Situations, the reference state of motion may be taken as one in which the 
spacecraft (solid + liquids) undergoes a uniform steady state rotation as a 
single rigid body about its designed spin axis under the influence of a con- 
stant thrust directed along that axis. With this approximation to the 
actual reference state of motion, and assuming further no velocities of 
liquids relative to the cavities, it follows that 


!- 

K . 

i 

f 



i 


w 

= mk , 

X S* 

0 , 

(4) 

x (U'to) 

a 

= ak , 

F - 

Fk 

, 

W X (I*W 






N 

U) 

- 0 , 

r = 
c 

= r 

c 

k 

• E 






n=l 


w] , 


-M-jO , 


( 2 . 1 ) 
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Ev 


, x [ (K-I) • .=] , h - K - - (M+Mj) ir^TE - r c "r c 


n-1 


K - I 


>[i* 0 i 2 * - Vc] ’ 


£ / [ i5 ni 2 E - V„] V“n ’ W J » 

n"l 


v„ * 0 , f, ■> 0 , r, m 2 ** , p * tu> ( f ~z ) , 

*' « ” * 1 n M * ** *<% * * 


n 


n n n rf / Q n » t n> 


V ■ u + .< x R , a * u + 7> x (~xR ) , 

n n o v n * 


(n“l »2j.,,,N)} 


( 2 , 2 ) 
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f - 0, 
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1 1 

tX 


| C t, + I"' 2 [( I l„ + >'n) 2+ ( L 2„%) 2 ] j + '-3n " 
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\ gn 



if 

‘W I 

4, 

(I '2n +y 

" g n 

*V ” 
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8 n n 

nr 


“ 8 n 

(x > 

n 

4* 
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>? 

H 
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4 



'n 


Q f o- 
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) + (L_ +y ) 

i » n* 


1 1 1/? - 
2n V J j 


(n*l , 2, . . . ,N) ; 


(2.3) 
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(2,4) 


Ma - P + 
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n*l 


n 


N 

I *- “ + ~ x (!•") - T + V 7 


n=l 


a V p . 

n f> n *n 
n 


V • V « V • v - 0, \ Q k i , 

m n n m 9 / **•» m ' 


n n n n 


V x a * 0 , 

n n 


v • v 
n n 


n n 


Q 4 o 
H n w 


f cos v ,k )» (f *0), Q * o r , (z *f ) 
n x n' n ' n ( F n n 

t C n 


p * 0 , Q r o„ , (z = f ) , 

r n F n n 

n 

F * I v p da * J v P do + / 
n n n w„ J n n w„ JL 


(n«l,2, , . , ,N) 


w 


do + 

" % " '' W " "°F 

n n 


v p do„ 
n r n F 


f - d f , 

I v p do * / V p 
/ n'n n 4 n^n 


dt 

n 


n 


C = / (R xv ) p do = / ( R xv 

n y„ n n n w n •'o n n 


w 


) P do„ 

n w 


w 


/ ( R xv 


J P n do F 


n •'of 


/ (R xv ) p do 
n n n n 


i* n 

[ - ( 

I (R xV P ) dx , / f r - 

J n n n n J n cos ^ v , 


n 


0=0 + 0 „ . 
n w F 
n n 


. ___ - ,k ) do„ = 0, 

Ot, t ti n F 

F n 

n 
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The constants e n appearing In (2.32) Cfln b® determined from the known liquid 
volume in each cavity. 

The foregoing completely define the reference state of motion of space- 
craft and liquids. 


PERTURBATION EQUATIONS 

In the reference state of motion defined above, the spacecraft with 
partly filled liquid cavities can be associated with a certain solid vehicle, 
called transformed consisting of the given vehicle and frozen liquids with 
free boundaries (2.3). I.et us find the small deviant motions of the trans- 
formed spacecraft in passing from a position corresponding to steady state 
rotation under constant thrust to a neighboring perturbed position. One can 
conceive of this transition as proceeding in two stages: 1) displacement of 

the whole spacecraft, as a single rigid body, to a perturbed position; 

2) deformation of the forms z n * f n (n»l ,2, . . . ,N) of the liquids (by cover- 
ing the free boundaries with layers of zero volumes) into forms 
zn * fn + *n (n*l ,2, , . . ,N) , where 


f * df , (n*l ,2, . . . ,N) (2.5) 

n n 


are small deviations relative to the f n , and i,\ denote differences in the 
forms of the instantaneous and reference volumes. 

Regarding the nth cavity, let,^ n stand for the region of the x n y n ~ 
plane bounded by the projection of a closed curve representing the locus of 
the Intersection points of the reference free surface with cavity wall(s) 
•i Wn , Then, in the first approximation Ref. 4, 



where (••■) can be a scalar, vector or dyadic point function. In particular, 
we observe 



f' dx dv * 0 
n n * n 


(2.7) 
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owing Co the incompressibility of Che liquid. Consider the disturbed form 
of F n , namely, 


F + F‘ 
n n 


/ 


/V p+V pM di 
\ n n irni 


( 2 . 8 ) 


n 


t +t ' 
n n 


I V p dT + / V p ' dt + If dx dy / 

J n r n n J nn n J J n n J 

J n l n *^n f 


rr J, 

ff f 

JJ dx n dy n J ''*K iz n ' 


f_+ f ' 
n 


V p dsi 
n r n n 


n 


f 

,/i n n 

In view of (2.6), where F , p , i are reference motion values, and 
' n *n ’ n 


n 


( F' » dF 
In r 


P* * d P„ 
n n 


(2.9) 


are small deviations. 

Subtracting out the reference value of F , i.e., 


/ 


V p 
n r n 


d P 


n 


n 
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from (2*8), it follows 



which, to the first order of small quantities, becomes 



But, from (2.2^) and (2,4^), 

V p = - p [ « + lit x 0»s x R )] 

n r n n n 


Substituting (2.12) into (2.11) and performing the indicated 
one obtains, to the same degree of approximation. 


F' 

n 


V p ’ dr - a) x (m X 

n n n 


f - - 

/ V p' dx - ai X [a) x 
J n r n n 

x 

n 




p dx dy ] 
n n 


^n 



dx dy ] 
n ^n J 




in view of (2.7). 


( 2 . 10 ) 


( 2 . 11 ) 

( 2 . 12 ) 

Integration, 

(2.13) 
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Using similar arguments, it can easily be shown that 


/ 


t: - / (R x V V)dT 


n n‘n n 


/ / ( R X u) p f ' dx dy 
JJn n n n^n 


Ik 


X [S X Q X R )] f ' P dx- dy , 
1 n J ) n n n •'n 


(2.14) 


c M+M 


l 


1 

M+M 


E ll R F dx dy 
JJ nn n n'n 


1 


n-1 


t If 


r f* p dx dy , 
n n n n n 


n=l 'j# 


(2.15) 


K' - 


i If 


[ | R ! \ *• R R ] f ' p dx dy 
1 1 n 1 n n J n n n 3 n 


o=l 


(2.16) 


In arriving at these results terms involving products of deviant quantities 
and/or their derivatives as well as local positional changes of liquid 
particles were neglected. 

Working to the same degree, of approximation. 


V 

n 


u' + <7»' * R + v' , 

n n 


V ' = u ' + m ' x R + ui x v ' + v ' , 

n n . n n fc 


V ' + m ' x V + m X V ' , 

n n n 


a ' + V) ' x R +• w 1 X (it' X R ) + >it X (us ' X R ) 
n n n 


+ 2u> X v' + v' 

n n t 


(2.17) 
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where 


a' * u' + m' x u + n' X u' 


(2.18) 


The perturbational equations that govern the small liquid motions in 
the N cavities are obtained in a similar fashion. They are 


a' « -~i v P' , 
n p n n 
n 


V . v' = V • v’ * 0 
n n n n 


V x a' - 0 , 

n n 


v • v ' 
n n 


Q n t. t , 
n n 


Q n f ' w » 
n 


(2.19) 


>(n=l , 2, • • *,N). 


f ' (v • k ) 
n n n 


Q n ' °F ’ W J 
n 


p' - P u f = 0 , Q r o , (z *f ) , 

r n n n ai Inn 


/ 

V. 11 F 


f coa(v , k ) da T , = 
n L n n I* 

t n 


//■ 

*^n 


f dx dy =0 
n. n 'n 
t 


To the same degree of approximation, it can be shown that 


, Q *■ a 
n w 


n 


( 2 . 20 ) 


vi • v' 
n n. 


(n=l,2, • • ’,N) , 


f (v • k ) , Q i: ( 1 „ , (z =f ) , 

n. . n n n F n n 

tt n 


and 


v • v' 
n n 


ttt 


n 


Q z a 

w 


( 2 . 21 ) 

(n=l,2,**-,N). 


fl (v„ • it,) , Q n g Op , (z„*fj , 


tttt 


n n 


n n 


The disturbed form of the translational equation of motion of the 
spacecraft is given by 

N 

M (a + 5 ') = ? + ?' + (? n + F^) (2.22) 

n=l 

where, as before, u, F, F are reference motion values given by (2.1) and 
(2.4g), and n 

da 

d¥ (2.23) 

d? 
n 

are small deviations therefrom. In the reference state of motion 



Ha = F + 



(2.24) 


in view of (2.4^). Subtracting (2.2^) from (2.22), one obtains 


Ma' = F' + 



(2.25) 
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for the perturbational equation of translational motion of the vehicle. 
Now, from (2.13), 




N N - 

£*:- £ / 

n*l n«l 


V p' di - «' x x 

n* n n 


E //' 


R f’p dx dy ] 
n n n n n 




N , 

E / V p' dt - (M+MJ m X ("> X v') , 

J n 1 n n 1 c 


(2.26) 


n«l i 


in view of (2.15). In addition, from (2 . 19^) and (2.17^), 

V P ' * - o ! n ' + m * X R + ' X (m X R ) 

n*n n | n n 

+ ut X (u\ 1 X R ) + 2«»» X V 1 + V ’ I , 
n 51 n t ) 


so that 


N r 

£ / 

n«l x 


V p' dr 
n‘ n n 


- M, n 1 


(M+M. ) m ' X r 


- (M+ML) m' X (io X £ ) - (M.+M) m X (ui' x r ) 
X c x c 



/ (v 1 + 2« X v') 

/ n * n 

J x C 

n 


n 


dr 


n 


(2.27) 


= - M «' - (M fM) m' X £ - (1L+M) m X (Hi' X £ ) 

IX c x c 


- £ 

n =i 


/ 


(V 

n. 


+ 


2 u* X v') 
n 


p dt 
n n 


> 




2-10 



where 



£ 




r it 
c 


W “ 03 


k 


Owing to (2.26) and (2.27), equation (2.25) becomes 


(M+M.) [ex' + u' X r + w X (u' X r ) + u x (u X r')] 
1 c c . c 


N 


£ j % 

n-l J x 

n 


+ 2w x v') p dt = F’ 
n n n 


Equation (2.29) can also be obtained by perturbing (1.28j) directly. 

The disturbed form of the rotational equation of motion of the 
spacecraft can be written 

N 

I ’ (w + U') + Q + U') X [I • (u + m')] = T + T' + ^ (T n +T^ 

n=l 


which, upon discarding nonlinearities involving m' and subtracting out 
reference equation ( 2 . 4 ^ ) » yields 


(2.28) 


(2.29) 


(2.30) 


for Che perturbations! motion. But, from (2.14) 


N 

E 

n-1 


N f 

E / 


(R XV p') di 
n tv n n 


n-1 i 
N 


n 


T, If 

n-1 jj> 

,n n 


t lf\'. 


(R X «) p f' dx dy 
n n n n ' n 


(2.32) 


x[mX(mxf)] r,' n dx n dy n 


n-1 .*? 


n 


E f ®n * VA 1 d \, ’ (M+ V ^ * “> 

n-1 ■'i 


(K 


/ ! 


uO 


i 


in view of (2.15), (2.16) and the identity 


i If 

n-1 >/> 

n 


R X |“ X (“ X R )] I f ' 0 dx dy » K' . 
n n n n n mi 


(2.33) 


Moreover , 


N r 

E / 


(2.34) 


(R x V p') it 
n n n n 


n-1 t 


n 


(M+M l ) (r x - (K-I) 


X [ (K-I) • .»] - x t (K-I) * 
N 


s/ 

n-1 t 


i R X (v ' 4- 2 X v ' ) ) P d t , 

| n n t n | n n 
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in light of the identity 


E 



R* n X [“' X (7x R n )] 


dx 


n 


♦ E 



R X 
n 


W x O' x 


R )] 
n 1 


n 


dt 


- X ( (K-I) • 7 ] + 7 X [(K-T,) • 7 '] 


( 2 . 35 ) 


Because of (2.32) and (2.34), equation (2.31) assumes the form 


K • 7 ' + 7 ' x (k . 7 ) + 7 x (k • 7 ') + 7 x (K' • 7 ) 

+ (M+M ) O' X a) + (M+M.. ) (r X a') 

Jl C X c 

N 


(2,36) 


£/ 

n*l t 


[R x (v' + 27 x v')] 0 dt « T' . 

n n„ n n n 

t 


Equation (2.36) can also he deduced by perturbing (1.28^) directly. 

The foregoing results are predicated on the assumption that v ,he 
steady state spin and all rotational and translational perturbations are 
instantaneously communicated to the liquids inside the cavities. As noted 
in Ref, 2, this assumption is at variance with the mechanisms of vortlelty 
transfer from the boundaries of the tanks to the contained liquids. This 
can best be understood by the following. 


Suppose that the cavities are motionless for t * t Q and that the 
rotation of the vehicle about its spin axis attains its steady state value 
by accelerating smoothly so as not to cause any ripples on the free 
boundaries. Then, under these circumstances, vorticity is, for the most 
part, transmitted by coriolis and centrifugal effects in thin Ekman layers 
(Ref. 5). On the other hand, these layers do not develop for periodic 
motions at spacecraft control frequencies. Consequently, even though 
steady state spin equilibria may be attained throughout the viscous liquids 
after some transient time, small periodic rotations about local axes of the 
tanks can never be transmitted to the cores of the liquids and should, 
therefore, be neglected. 
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I 

t 

i 

The above arguments require that the terms >.i' * r n , x r n in (2.17) 
be discarded. With this injunction, the perturbational equations that govern 
the motions of spacecraft and liquids become 


N 


Ma* - P* + 2 ? n ’ 

n*l 


N 


I • u»' + .o' X (I • w) + t.1 X (I • ;.' f ) * T 1 + 2^ 

n*l 

1 


a ’ * V p 1 , 

n p n n H n 


•\ 


v • V' = V • v' ■ 0 , Q i: T , 

n n n n n n 


V a' * 0 , 

n n 


v • V = 1 
n n 


, Q v a 
^n w 


| f' (v • k ) , Q v. o„ , (z *f ) , 

t n ' n n ai F.. n n 


p 1 - no f * 0 , Q c , (z =f ) , 

1 n n n ni F n n 


n 


I f ' cos (v , k ) da„ = II f * dx dy 

J n t n n F n JJ n t n n 


t# r 


(2.37) 


F' = I V p' di - rn x [w x /./ R f' p dx dy ] 
n J n n n J J n n n n •'n 


>(7=1.2, • — ,N), 




n 


/ v p diT, 
n n w 

n 


ui x [w X 


II R f P dx dy ] , 

J J n n n n J n 


w 




n 


n 
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SECTION 3 


PERTURBATIONAL LIQUID MOTIONS REDUCED TO SOLUTIONS OP N 
INHOMOGENEOUS BOUNDARY VALUE PROBLEMS 


PERTURBATIONAL EQUATIONS DESCRIBING LIQUID MOTIONS 

In light of (2.37) and (2,38), the perturbatlonal aquations that govern 
the motions of liquids Inside the N cavities can be written 


r 


v* + 23 X v ' + V i' - o , 
n„ n n n 

t 


\ • \ ■ °< 


(’ + 2a> X 5' - 27 (ui-v’) - 0 , 

v n_ 'n n n 

t 

<f>' * p'/p„ + B n • R rt » 
n n n n n 


B-a'+S'XL + m' X (~XL ) + w X (3 'xl ) 
n n ' n 


w * wk * wk n , constant) 


v • v 
n n 


» - 


.Q et 
,x n n 


» Q £ a,, , 

n w ’ 


f' (v *k ) , Q e o_ , (z - f ) 
n ' n n F n n 

t n 


p' - ap f' = 0 , Q e cr_, , (z * f ) , 

K n ^n n H n F n n 


n 


j \ <V E o> d \ ‘ [fa dx n% 

\ K 


- 0 , 


(3.1) 


V (n-1,2, . . . ,N) 
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In addition* 


v • v 
n n. 


0 , 


f'„ (v • It ) . 

n tt n " 


Q to , 
n w 

n 


Q e o„ , (a * f ) , 
F n n 
n 


v v 
n n 


ttt 


0 , 


tttt 


Q » o , 
n w 

u w tx 


E „ ... < 7 n • V - V V 


n 


v„ • - v„ • v;^ - 0 , Q n T n . 

[ C n < 7 „ ' V 6 \ ■ JJ E n dx n djr n 


0, 


n 


I, 


l 


dx - US X [7,1 X II R £ 0 dx dy 1 

n if n n n n n 


v p do - » x [u * R f ' p„ (v * k ) do ] , 

n r TtW * nnnnii l 1 


n n ' n n' F 


w_ 


n 


'Op 


n 


n 


f <\ * Vn> 
J T 


dx 


{R x [a + u x (u x R )]} 
n n 


R 


n 


f p dx dy * | (R x v ) p af 

n n n n f n n 

' V 

J W 

n 


/ tR n* 


[a 


+ "J x ( M x R n )]> f n P n (\> n * k n ) dOp , 

n 

(v„ • k„) do^ - / p n dx n dy n - 0, 


p \ v Iv ) uu 

*\n n n F_ 


vi 




n 
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DERIVATION OF BOUNDARY VALUE PROBLEMS 

The determination of perturb a tional liquid motions Inside the N cavities 
can be made to depend on the solutions of N boundary value problems with 
acceleration potentials as dependent variables. This can be shown In the 
following manner. 


Consider the nth cavity. Take the divergence of (3.1^, 


1 • /v» +2uXv' +v n«0 , 

n I n^. n n T n J 


obtaining 




(3.3) 


in view of 


V • v' 
n n. 


0 , V • V <fr' 

n n n 


V V , V n • (25 v*) 

n n n n 


2w • V 


Differentiate (3.3) partially with respect to time, 


vV = 25 • ?' 

n t n t 


(3.4) 


Take the -fcalar product of to and equation (3 5 

-i. 

= 0 , 


2w 


L t 


+ 2(0 X 5' - 2V (wv') 

n n n 


obtaining 


2m • £' = 4m • V (m*v') , 

n n n 


(3.5) 


in light of 


2m • (2m xg =0 
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Substitute (3.5) into (3.4) and differentiate the resulting expression 
partially with respect to time, 



(3.6) 


Take the scalar product of m and equation (3.1j) 

m • I v ' + 2w X v ' + V (t ' I ■ 0 , 

n„ n n r n ’ 


obtaining 


v 


- i,* 


V ip’ 
n n 


(3.7) 


owing to 


tii • (2t«j X v') =0 . 
n 


Substitute (3.7) into (3.6), 


vV + 4u> 2 « 0 , Q t* r 

n n. . n „ hi n 

tt z z 

n n 


In view of 


m • V (iii • V d> f ) = 

n n n n 


2 2 
n n 


(3.8) 


Expression (3.8) is the desired partial differential equation for 

To establish boundary conditions for first differentiate (3.1^) 
partially with respect to time. 


v ' 


+ 2ui x 


tt 


n. 


V if)' 
n iv 


(3.9) 
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Now, from (3.1^) 


- V if)' - 2u X v' , 


n. n'n 

t 


(3.10) 


Substitute (3.10) into 3.9), 


\ t + ^ - '"‘'X + \\ + V „ X ^ - 0 ' 


(3.11) 


in view of 


V x (2“i(j.') « - 2to x V if)’ . 
n n n n 


Differentiate (3.11) partially with respect to time, 


v' + 4 io 2 v' - 4to 2, 
"ttt n t 


ic (k -v W V <}> ’ + V x / 2to<j> ' \ = 0 . (3.12) 

n \ n nj n^ a ( ’nj ' 


Take the scalar product of and equation (3.1^), 


v • k 
n t n 


4 > ' 
T n 


n 


(3.13) 


because 


k • (to X v') 0 , 

n n 


k • V <} ' = <b ' 

n n n T n 

z 

n 


Substitute (3.13) into (3.12) and rearrange. 


V <)>' + V 

n n tt 


11 l N 


+ 4to d> ' k = 
r n n 
z 

n 


/v' +4<o 2 v' \ 

\ n ttt n t/ 


(3.14) 


Project (3.14) along the outward directed normals to a and , 


n n 


5 • V A 1 + V X /2w 4> * \ 

n n’n tt u ^ nj 


+ 4m $ k 
^n n 
z 

n 


“ 0 , Q i cr 
n w 


(3.15) 


*-(v • k ) IV + 4m 2 £* \ , Q t o » (z “ f ) » 

n n \ Vtt n tc) n F n n n 

In view of (3.2^) and (3.22). To complete the specification of the boundary 
conditions, solve (3 . 14 ) for p^ 


p 1 * p (<Ji ’ - 0 » R ) 

*n n n n n 


and substitute Into condition (3.1g), obtaining 


<ft* - 
T n 


of' 

n 


R_ 


Q n W F 


(3.16) 


The free surface displacement can be eliminated from the problem as follows. 
Solve (3.16) for F and substitute the resulting expression into (3.15), 


v • V <j>’ + V 

n n r n n 



4<»w(> ' 
n 


n 



= 0 


Q 


n 



(3.17) 



L' 

\ n tttt 


4m^tj) 


R 


tt , 


Ik h 

\ tttt 


+4m 2 (i 


tt 


Q cct 

ui F 


(z«f ) 
n n 


where 


f ’ 
n 


— 4>* - - 0 
a T n a n 


R » Qn e S 
n n r 

n 


g 


n 



(Equation (2.3^)) , 


Q c 


F 


(3.18) 


(3.19) 


3-6 


Take the divergence of (3.14) and make use of (3.2^), 


V • 

n 

VJI + X 

n n n 

tt 

|2m 

\ V 

\+ 4w^' k 
) T n n 

' z 




n 


V 2 d>' 

n r n 


tt 


4m ^<|) ' 
n 


z z 
n n 


Q e 

^ n 


x 

n 


> 


(3.20) 


in agreement with (3.8). 

The inhomogeneous boundary value problem defined by (3.8) and (3.17) 
governs the perturbational liquid motions inside the nth tank. The effects 
of centrifugal and corlolls accelerations together with vortlcity are 
implicitly taken into consideration in tie formulation. 


Observe that the excitation vector Bn is a function of the perturba- 
tional translational and rotational accelerations and rotational velocities 
of the spacecraft. They, in turn, are related to the motions of the liquids 
in the remaining N-l tanks via the perturbational equations of translational 
and rotational motion of the spacecraft. 


The force and moment resulting from the action of the perturbational 
liquid motions on the tank wall can be written 


F’= vfp da - / 

n J n T n n n J 


a; x (m x R ) 
n 


<£’ p da^ - 8 
T n n F n 
n 


/- 


w x (w X R ) _ _ 

. — ^ • R ) p da„ 

g ' n n n F 


n 


n 


Ip dt 

J n 11 


n 


(3.21) 


T' = 
n 


f _ f I R X [a+wx(wXR )]i 

/ (R X V ) I' p do - / 

I s n n Y n f n n J g n T n 


K P n d<, F 


n 


n 


- /«„ X 


n 


n 


B) P n dx + 

n n n 


/ 


{r r X [a+wX(wXR n )]| 


(3.22) 


J n 


(B *R ) p da„ , 
' n n n F 


n 
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ST’ 


In light of (3 . I 4 ) ( 3 . 25 ) (3.2g), ( 3 . 27 ), (3.18), (3.19) and the divergence 
theorem. 


VARIATIONAL FORMULATION OF ASSOCIATED HOMOGENEOUS PROBLEM 

Consider the associated homogeneous boundary value problem pertinent to 
the nth tank. 


V + V 
n r n. . 
tt 




+ 4u * ic 
n n 
z 

n 


V <f>* 
n n 


+ V 


tt 


X /2wij> ' \ + ' 

n [ \l 


n n 

7 

n 


0 , Q n e T n , 


0 , Q e. 0 , 

^n w 

n 


(3.23) 


(v -it ) 


^ (♦;«« + ) • ^ c \ • (z - ■ y 


It is fairly straightforward to show that the first variation of the 

functional (^1 = 0 ) 

n 



n 



(3.24) 



subject to 


*n 



t 


0 , 


gives the boundary value problem (3.23). 
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The functional (3.24) together with the methods of Rayleigh and Rltz 
make It possible to obtain approximate solutions of (3.23) in terns of 
eigenfunctions and eigenvalues of simpler, but related, boundary value 
problems. 

For Instance, the eigenfunctions - 4 = and eigenvalues K satisfying 
the boundary value problem n i n i 


v \ - 0 • e 'n ' 


V • V 1 |' 

n n n , 


Q r 0 , 

n W n 


(3.25) 


K 

Jjl 

n i ’ 


v) i, o,, , (z ■ f ) 

<n F ' n n 
n 


and orthogonality conditions 



where 


are available for a spherical container, Ref . 1. Assume that <j^ can be expanded 
in n series 



where, q n ,(t) are to be determined. Substituting (3.28) into (3.24) and 
performing the first variation of the resulting expression (51^ « 0) 
subject to 


5 q 


n. 



t 


m 0 


t 


one obtains the system of ordinary differential equations 



for the determination of the. q . Other examples can be adduced. 

”i 

If the excitation vectors are regarded as preassigned functions of 
time, the inhomogeneous boundary value problem defined by (3.8) and (3.17) 
can be deduced from the first variation of the functional 
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CHARACTERIZATION OF THE BOUNDARY VALUE PROBLEMS 
For a single frequency excitation 

(3.31) 

the inhomogeneous boundary value problem spelled out by (3.8) and (3,17) 
becomes 
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(3.32) 


n 


3i|j / , 2 \ Dip „ 

n t . n t . 1 , 4m 1 n . 2m , 

~ — x + T — J +11 r- I- — k - x -T- V x k i 

3x n 3y J n I „2 I 3z n S n n n 

n / n l 52 / n 


= 0 , 


Q c a 
n w 


(3.33) 


n 


si/i-iai)* « 2 


4 w* 


g. 


n 


52 


n 8. 


(c *R ) , Q 
n n n F 


n 


3-11 


If 2.i/:: <1, the partial differential equation (3.32) can be brought to an 
elliptic form by the transformation 
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(3.34) 


with boundary conditions 
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whore t n , o Wn , n F arc* the stretched values of T n » o Wo » UF n * On the other 
hand, if 2'V./u >1, the partial differential equation (3.32) can be brought 
to an hyperbolic form by the transformation 
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with boundary conditions 



For low spin rates the liquid motions are governed by an elliptic 
partial differential equation with mixed boundary conditions. As the spin 
rate is decreased indefinitely the governing equations degenerate to the 
classical lateral slosh equations. This strongly suggests that the custom 
of representing sloshing motions in spinning tanks by simple mechanical 
pendulums is valid only for 0 < 2 m /'l < 1 . 

For ratios of 2,/- in excess of unity, the sloshing motions are governed 
by an hyperbolic partial differential equation with mixed boundary condi- 
tions. The nature of the hyperbolic differential system together with the 
form displayed in (3.19) point out that the representation of sloshing 
motions in spinning tanks by mechanical pendulums is, in general incorrect. 

Yet, for most practical applications 2h\/u >1. Consequently, the use 
of mechanical pendulums to reproduce sloshing motions in spinning tanks is 
suspect and should be avoided, unless severe restrictions can be tolerated. 

The nature of the solutions of the above boundary value problems in 
the case of a completely filled spherical tank is discussed in Ref. 5. 
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SECTION 4 


SUMMARY OF PERTURBATIQNAL EQUATIONS 


The perturbations! equations that describe the small motions of a space- 
craft with N partly filled liquid propellant tanks relative to a reference 
state of motion wherein the vehicle undergoes steady state rotations, as a 
single rigid body, about its designed spin axis in a constant acceleration 
field are collected below. 
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and the reference state parameters are spelled out in ( 2 . 1 , 2 , 3, 4), 


CONCLUDING REMARKS 

The perturbational liquid motions in the nth cavity are governed by the 
boundary value problem displayed in (4 . 13 ) and ( 4 .I 4 ). Further meaningful 
progress is not possible until a method is devised for solving this differ- 
ential system. 

For a single frequency excitation and low spin rates, the. problem reduces 
to an elliptic partial differential equation with mixed boundary conditions. 
This differential system can be solved numerically in a fairly straight- 
forward manner. However, this case does not seem too important in most 
practical applications. 

For a single frequency excitation and high spin rates, the problem 
reduces to an hyperbolic partial differential equation with mixed boundary 
conditions. Numerically, the solution process of this differential system is 
formidable. However, a numerical approach using Green's function appears 
promising. This would, of course, require extensive analyses to develop the 
computational algorithms. 
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Another approach to the general problem is to employ a variational 
formulation in conjunction with the methods of Rayleigh and Rltz. Basi- 
cally, the technique depends on being able to represent </>n In the form 


fr' 

n 



q 


(t), 


where the %£ ai*o eigenfunctions of a related, but simpler, boundary value 
problem, and the q^ are generalized coordinates to be determined from the 
solution of an infinite set of fourth order ordinary differential equations. 
This set of equations is obtained by minimizing the variational Integral 
with respect to the Independent coordinates q n i* Such a method was used 
to arrive at system (3.29). In light of the discontinuous nature of the 
problem, this approach should be used with caution. 
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